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, ^ Abstract. We present the results of a massive variability search based on a photometric survey of a six square 

^ . degree region along the Galactic plane at {I — 305°, b = —0.8°) and {I = 330°, b = —2.5°). This survey was 

performed in the framework of the EROS II (Experience de Recherche d'Objets Sombres) microlensing program. 
The variable stars were found among 1,913,576 stars that were monitored between April and June 1998 in two 
passbands, with an average of 60 measurements. A new period-search technique is proposed which makes use of a 
statistical variable that characterizes the overall regularity of the flux versus phase diagram. This method is well 
suited when the photometric data are unevenly distributed in time, as is our case. 1,362 objects whose luminosity 
varies were selected. Among them we identified 9 Cepheids, 19 RR Lyra, 34 Miras, 176 eclipsing binaries and 266 
Semi-Regular stars. Most of them are newly identified objects. The cross-identification with known catalogues 
has been performed. The mean distance of the RR Lyrffi is estimated to be ~ 4.9 ± 0.3 kpc undergoing an 
average absorption of ~ 3.4 ±0.2 magnitudes. This distance is in good agreement with the one of disc stars which 
contribute to the microlensing source star population. Our catalogue and light curves are available electronically 
from the CDS, Strasbourg and from our Web site". 
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1. Introduction 

In 1996, the EROS II collaboration started an observa- 
tion program towards the Galactic Spiral Arms (GSA) 
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dedicated to microlensing events. Since then, four regions 
of the Galactic plane located at large angles with respect 
to the Galactic Centre are being monitored to disentangle 
the disc, bar and halo contributions to the microlensing 
optical depth. Seven microlensing event candidates have 
already been published, based on three years (1996-98) of 
observations (Derue et al. 1999, 2001) [hereafter papers I 
and II] . The distance of the source stars used in these pa- 
pers to compute the expected optical depths was deduced 
from a detailed study of our colour-magnitude diagrams. It 
was thus found that the source star population is located 
~ 7 kpc away, undergoing a n interstellar ex tinction A{V) 
of about 3 magnitudes (see Mansoux (1997 ) for more de- 
tails). This distance estimate is in rough agr eement with 
the distance t o the spiral arms ob tained by Georgelin el 
al. (I9Q and Russeil et al. (1998 ), but its uncertainty is 
limiting further interpretation of our microlensing optical 
depth estimates. It was therefore desirable to seek more 
information on the distance distribution of the source stars 
- whether these stars belong to the disc or to the spiral 
arms - and on the reddening along our observation line 
of sights. This led us to perform a dedicated variable star 
search between April and June 1998, on a subset of our 
Galactic plane fields. The analysis was restricted to the 
brightest stars of this subset. 

Among the wide variety of variable stars, periodic ones 
are of particular interest. The properties of Cepheids make 
them well suited to trace the Galactic spiral arms. Their 
reddening is measurable as well as their distance via the 
period-luminosity (PL) relation. RR Lyrae stars are old 
stars, well suited to trace the disc population. One can 
infer their mean dereddened magnitude a nd their abso- 
lute magnitude ( [Gould and Popowski 1998 ). The infrared 
PL(K) relation for Miras and Semi- Regular variable stars 
can be calibrated using a comparison of DENIS and 
EROS LMC giant stars (|Cioni et al. 200l|) . Finally de- 
tached eclipsing binaries also offer the opportunity to mea- 



sure their stellar parameters and their distance (Paczyiiski 
19960 

This paper presents the results of this particular cam- 
paign which led to a catalogue containing a large number 
of new variable objects in the Galactic plane. Sect. gives 
the basic features of the observational setup. Sect. |^ gives 
details on a new algorithm used to search for periodic vari- 
ations of the luminosity. Sect. ^ describes the catalogue 
and the cross-identification process. In Sect. I we use the 
selected RR Lyras to estimate the mean reddening of our 
fields and we give the distance distribution of these stars. 

2. Experimental setup and observations 

The MARLY telescope and its two cameras, the way we 
carry out our observations, as well as our data reduction 
sequence are described in Paper I and references therein. 
The two EROS passbands are non standard. The so- 
called EROS-red passband Re is centred on A = 762 nm, 
close to Ic Cousins, with a full width at half maximum 
AA ~ 85 nm, and EROS-visible passband Ve is centred 



on A = 600 nm, close to Vj Johnson, with AA ~ 78 nm. 
The EROS II colour magnitude system is defined as fol- 
lows : a zero colour star with Ve — Re = (a main 
sequence AO star) will have its Re magnitude numeri- 
cally equal to its Cousins Rc magnitude and its Ve mag- 
nitude numerically equal to its Johnson Vj magnitude. 
The colour transformation between the EROS II system 
{Ve,Re) and the standard Johnson-Cousins {VjJc) sys- 
tem is then : 



Ic — Re 
Vj ^ Ve- 



- 0.01 X {Ve - Re) 
0.39 X {Ve - Re) ■ 



(1) 



The colour coefficients are obtained from the study of 
our passbands based on Landolt standards and on one 
of the EROS II LMC fields observed simultaneously in 
BVjRcIc with the Danish 1.54 m (at ESO-La Silla) and 
with the MARLY. The zero points are established with 
t ertiary standar ds in BVRI taken with the Danish 1.54 
( Regnault 2000 ). We have cross-che cked our Re photom - 
etry with the Ic one of DENIS ( |Fouquc et al. 2000| ). 
Furthermore, using Eqs. (1), the mean magnitudes of the 
LMC red-giant clump stars agree within 0.1 magnitude 
with determinations made by Harris & Zaritsky (1999) 
and Udalski et al. (1998). We thus estimate that the pre- 
cision of the zero points of the MARLY calibration is 
~ O.l'"'*^, 



Table 1. Characteristics of the 6 one square degree 
fields monitored for this study. The equatorial (J2000) 
and galactic coordinates of the field centres, the number of 
photometric measurements per light curve per colour N„i 
and the number of analysed light curves Nobs (in millions) 
are given for each field. The observation duration was ~ 
100 days. 



Field 


Q 


5 


6° 


r 




Nobs 


Norma (7 Nor) 


59 


1.30 


gn450 


16:09:45 


-53:07:03 


-1.17 


330.49 


63 


0.28 


gn453 


16:22:28 


-52:06:20 


-1.69 


332.24 


52 


0.36 


gn455 


16:26:52 


-52:21:02 


-2.35 


332.54 


57 


0.30 


gn459 


16:15:51 


-54:48:45 


-2.86 


329.82 


61 


0.36 


Musca {e Mus) 


60 


0.61 


tm550 


13:27:04 


-63:02:18 


-0.47 


306.98 


61 


0.27 


tm551 


13:31:18 


-63:34:41 


-1.07 


307.37 


60 


0.34 


Total 


1.91 



Among the 29 fields of the EROS GSA microlensing 
program, the six fields considered here were monitored 
about once per night between April and June 1998. They 
represent 2 square degrees towards 9 Mus and 4 towards 
7 Nor (named after the closest bright star). Table |l| gives 
their coordinates, the number of images iV,„ taken in each 
direction and the number of analysed light curves Nobs- 
To avoid CCD saturation by the brightest stars {Ic ~ 9), 
by Cepheids or Miras in particular, we have reduced the 
exposure time to 15 s instead of the 120 s used in the 
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and Ve- The photometric accuracy is ~15% at i?£;~18, 
and about 2% for the brightest stars. Finally, using the 



S 1C- 




4 116 118 120 122 124 126 128 130 132 134 

Week number since 01/01/96 

Fig. 1. Average time sampling for the 6 fields monitored 
(2 towards 7 Nor, upper panel, and 4 towards 6 Mus, lower 
panel), in number of measurements per week and per field. 



microlensing survey. As a consequence the catalogue is 
incomplete as far as faint stars are concerned, but it could 
be updated later by using the total set of available GSA 
images. Fig. 1 shows the average time sampling. Three 
gaps can be seen in our data : the first two (around weeks 
119 and 123) were due to bad weather conditions while the 
third one (around week 127) corresponds to the annual 
maintenance of our setup. 

3. Search for periodic stars 

3.1. Reconstruction of the light curves 

Since the EROS photometry is described in detail in 
Ansari (1996 ) only the main features of the PEIDAH — h 
package are summarised below. For each field, a tem- 
plate image is first constructed using one exposure of very 
good quality. A reference star catalogue is set up with 
this template using the CORRFIND star finding algo- 
rithm ( Palanque-Delabrouille et al. 1998). For each sub- 
sequent image, after geometrical alignment with the tem- 
plate, each identified star is fitted together with its neigh- 
bours, using a PSF determined on bright isolated stars 
and imposing the position from the reference catalogue. 
A relative photometric alignment is then performed, as- 
suming that most stars do not vary. Photometric errors 
are computed for each measurement, assuming again that 
most stars are stable, and parametrised as a function of 
star brightness and image sequence number. Fig. 2 shows 
the mean point-to-point relative dispersion of the mea- 
sured fluxes along the light curves as a function oi Re 




Fig. 2. Average value of the relative frame to frame dis- 
persion of the luminosity measurements versus Re (upper 
panel) and Ve (lower panel), for stars with at least 30 
reliable measurements in each colour. This dispersion is 
taken as an estimator of the mean photometric precision. 
The superimposed hatched histograms show the EROS 
magnitude distribution of monitored stars. 



PEIDA++ photometric package, we reconstruct the light 
curves of 1,913,576 stars. 



Table 2. Impact of each selection criterion on the data. 
For each cut, the number of remaining light curves is given. 



Cut 


Criterion 






Total analysed 


1,913,576 


1 


iVm > 30 


1,299,690 


2 


Re < 17 


330,089 


3 


Pre-filtering 


41,545 


4 


Period search 


2,553 


5 


Aliasing 


2,424 


6 


Visual inspection 


1,362 



3.2. Preselection 

Each one of the light curves is subjected to a series of 
selection criteria in order to isolate a small sub-sample 
on which we will apply the time consuming period search 
algorithms. These analysis cuts are briefly described here- 
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after (see Derue (1999) for more details) and their effect 
on the data is summarised in Table ^ : 

cut 1 : At least 30 measurements should be available in both 

passbands and the base flux must be positive; 
cut 2 : The search is restricted to stars whose magnitude is 

Re<^7 which corresponds to a photometric accuracy 

in Re better than ^^10%. 
cut 3 : A non specific pre-filter is applied which retains most 

variable stars. It selects light curves satisfying one or 

both of the following criteria : 

— the relative dispersion of the flux measurements is 
25% larger than the average one for the set of stars 
having the same magnitude; 

— the distribution of the deviations with respect to 
the base flux is incompatible with the one expected 
from a stable source with Gaussian errors during 
the observation period (Kolmogorov-Smirnov test). 

These cuts are tuned to select ~10% of the light curves. 
We have checked that this procedure allows one to 
retrieve the previously known Cepheids observed by 
EROS in the Magellanic Clouds. We also keep a ran- 
domly selected set of light curves (^2%) to produce 
unbiased colour-magnitude diagrams, for comparison 
purposes; 

At this stage a set of 41,545 light curves remains which is 
then subjected to a periodicity search. 

3.3. Light curve selection 

We use three independent methods to extract periodic 
light curves. The first two are classical methods al- 
ready described in the literature : method 1 is based on 



the Lomb-Scargle periodogram (see pcargle 1982| ) while 
method 2 mak es use of the One Way Analy sis of Variance 
algorithm (see 3chwarzenberg-Czerny 1996 ). Both provide 
the probability for false periodicity detection. In method 
1 one computes the Fourier power over a set of frequen- 
cies. It is therefore well adapted to identify sinusoidal light 
curves. It can be improved by incorporating higher har- 
monics in order to detect any kind of variability such as 
eclipsing binaries ( |Grison 1994 ; Orison et al. 1995 ). 

We developed a new method, the third one, in order to 
extract periodic light curves in a way which is insensitive 
to the particular shape of the variation. This method also 
provides a probability for false periodicity detection. 

It consists in searching for a frequency such that the 
corresponding phase diagram, i.e the series of fluxes Fi 
versus phases ipi in increasing order of ip, displays a reg- 
ular structure significantly less scattered than for other 
frequencies. Let Tabs be the observation duration (~ 100 
days in this analysis) . We span the frequency domain from 
T~f^l to (0.2 days)~^ with a constant step of (4 x Tots)^^- 
The total number of test-frequencies is thus Ntest ~ 2000. 
This sampling ensures that the total phase increment over 
Tabs is 7r/2 for two adjacent test-frequencies. 

For each value of the test-frequency we compute the 
corresponding phase diagram. We calculate a from the 
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Fig. 3. Distributions of (left panels) and variation of x 
with the value of the test-period in days (middle panels) 
obtained with the third method for a stable star and for 
a typical Cepheid candidate. The bold dots are pointing 
to the actual period of the star. The right panels show 
the light curve (in i?^;) obtained once the period has been 
folded in. 



weighted differences of Fi and the fluxes interpolated be- 
tween Fi^i and i^i+i : 



/V 9 

g fF, - (1 - R,) X F,^, -R,x F,+i ' - 
1=1 



(2) 



where R, = {(pi - ipt-i)/{(pi+i - <Pi-i) and is the 
number of measurements. The uncertainty ai takes into 
account the errors on the flux Fi and on the interpo- 
lated flux : af = ef -I- (1 — Ri)^ef_i + Rfefj^i. Expression 
(H) can be interpreted as the of the set of differences 
between the odd measurements with respect to the line 
joining even ones, added to the of the set of differ- 
ences between even measurements with respect to the line 
joining odd ones. If the star is measured in both colours, 
we add the obtained in each colour [Nm is then twice 
as large). For a given stable star with Gaussian errors, 
each phase diagram can be considered as a random real- 
isation of the light curve. When the test-frequency spans 
the search domain, the distribution of the x^ parameter 
defined by Eq. (2) is the one of the standard x^ with Nm 
degrees of freedom. Since Nm is large enough, this distri- 
bution is close to a Gaussian with average Nm and vari- 
ance 2Nm (see upper left panel of Fig. 3). For a periodic 
variable star, the x^ distribution displays a main cluster, 
when the test-frequency results in a phase diagram with 
non-correlated point to point variations, and a few lower 
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Fig. 4. Distribution of — \og{Prob) for method 1 (top), 2 
(middle) and 3 (bottom). The histograms show the distri- 
bution for all 41,545 stars on which the periodicity search 
is done (in white) and for the 2,553 selected light curves 
(in black). 



values when the test-frequency corresponds to a phase di- 
agram with a regular structure (see lower left and middle 
panels of Fig. 3). In practice, instead of using the param- 
eter defined by Eq. (0), we use the reduced variable : 



X 



<X^> 



1 



(3) 



where <x^> is the average of the realisations of for all 
test-frequencies. For a stable star, the distribution of this 
variable x is a Gaussian centred at zero, with unit vari- 
ance. If the errors are correctly determined < > /Nm 
is close to unity; if the errors are all systematically over- 
estimated (or underestimated), then including this term 
ensures a global renormalisation of the errors in Eq. (^, 
and the distribution of our reduced variable x is also a 
normal distribution. Let Xmin be the smallest value of 
X calculated among all test-frequencies for a given star. 
Under the hypothesis that the light curve is produced by 
a stable star, the probability to obtain at least one value 
X < Xrnin in a scrics of Ntest realisations is Prob{x < 

1 2 ,„ _ 

/2tt J-oc 



\Nt, 



1 



1 



If this 



probability is small, then 



Prob{x < Xrmn\Ntest) 



-erfc 



V2 



(4) 



If the light curve exhibits periodic variations, then there 
exist test-frequencies for which x is significantly smaller 
than typical values of this variable (see middle and right 



panels of Fig. 3), and the probability for false detection 
is then extremely small. Fig. 4 displays the probability 
distribution obtained with the three methods for the set 
of filtered light curve^ 

We apply the three algorithms which all give a prob- 
ability for no periodicity. A star is accepted only if se- 
lected by all three methods with the following thresh- 
olds : P(method 1) < 10'^^, P(method 2) < 10'^°, 
P(method 3) < 10^'^, tuned in order to allow one to re- 
trieve the previously known Cephcids observed by EROS 



in the Magellanic Clouds (see Deruc (1999) for more de- 
tails). This procedure {cut 4) selects a sample of 2,553 
stars. 

Ten times more stars would have been selected if we 
had used method 1 or method 2 only (with the same 
thresholds), most of them being spurious variables. The 



10 



10 10 

Period (days) 



Fig. 5. Period-colour diagram (P in days vs Vj-Ic) of the 
2,553 selected candidates before cut 5. 



third method would have added far less candidates if used 
alone, but still a factor 2 more. Combining independent 
methods has thus the advantage of considerably reducing 
the noise background (mostly due to aliases of one day or 
noisy measurements) while giving redundant information 
about the period. To obtain individual periods we perform 
Fourier fits with five harmonics : 



flux = fluxr 



,2tt. 



5 

-^a/ cos{^l{t ~ to) + (/);), 
1=1 



(5) 



The third method gives a probability which is always larger 



than 10" 



Indeed, the smallest value of x which can be ob- 



tained takes place when x ~ 0, then 



/iV™:< 10. 
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where P is the period, (pi the phase and a/ the amplitude. 
We define the amplitude ratios Rki = au/ai, and the phase 
differences (defined modulo 27r) ^ki = 4>k — k<j>i , with k > 
I. Objects with non-significant harmonic amplitudes (z.e 
with almost sinusoidal light curves) have i?2i ~ and 
their $21 is iU defined. 

The selection of periodic variable stars is complicated 
by aliases. Some of the stars with periods equal to a sim- 
ple fraction or a low multiple of one day may be badly 
phased because of the nightly sequence of measurements. 
These aliased periods are seen in Fig. 5 as vertical groups 
of dots at 2/3, 2 and 3 days. To eliminate them we de- 
mand {cut 5) that the fitted periods are not within ±1% 
of these values. One can also notice some vertical groups 
of points around 25 days which correspond to data gaps 
in our sample (see Fig. 1). Once these objects are re- 
moved, 2,424 stars remain. The flux values of the remain- 
ing stars are folded using each period obtained with the 
three methods. The resulting phase diagrams are visually 
inspected. Some of them display an obvious spurious peri- 
odic or quasi-periodic variability due to a low photometric 
quality. After this final visual selection {cut 6) the list of 
variable stars includes 1,362 candidates which exhibit un- 
ambiguous periodic variability. 

3.4. Type of variability 

The classification of the selected stars among different 
types of variability cannot be based on the position of the 
objects in the colour-magnitude diagram since the spread 
in distance of these stars entails a spread in magnitude 
and colour. It is desirable however to classify the vari- 
ous light curves according to some physical parameters. 
In the following we mainly use criteria based on the pe- 
riod P of the luminosity variations and on the amplitude 
ratio AV£;/Ai?E. For each selected type of variable star 
the phase diagram of a typical candidate is displayed in 
Fig. 6. Three groups are distinguished depending on their 
period : 



Table 3. Number of selected objects for each type of vari- 
ability. 



Period range 


Type 


Number of objects 


P > 60 d 




543 




LPV 


509 




Miras 


34 


60 > P > 30 d 




387 




Semi-Regular 


264 




miscellaneous 


123 


P < 30 d 




432 




pulsating 


60 




RRc 


14 




RRafc 


5 




classical-Cepheids 


6 




s-Cepheids 


3 




miscellaneous 


32 




non-pulsating 


372 




EA 


130 




EB 


35 




EW 


11 




miscellaneous 


196 



1*^' : sta: s with a period larger than 60 days (543 objects) : 
For these objects an entire period has not been ob- 
served. There is thus no warranty that these objects 
are periodic ones. 34 display a nearly linear light curve 
and are thus catalogued as Miras candidates. The other 
509 objects are catalogued as Long Period Variable 
stars (LPVs). 

2"^^ : stars with a period between 30 and 60 days (387 ob- 
jects) : 

264 Semi-Regular variable stars are selected by requir- 
ing : AVe/ARe > 1.2 to select pulsating stars (see be- 
low), Vj-Ic > 2.5 to discriminate from bluer variable 
stars and ARe < 1.0 to avoid possible Miras or LPVs 
wrongly phased. The long term stability of these stars 
is not known. Some of the reported periods may change 
from season to season, as a result of their semi-regular 
behaviour. The remaining 123 objects are catalogued 
as miscellaneous variable stars. 



3"^ : stars with a period smaller than 30 days (432 objects) : 
The colour change for a Cepheid in standard passbands 
is AVj/AIc >1.3 ( iMadore et al. 199l| ) which corre- 
sponds to AVe/ARe > 1.2 in the EROS system. Two 
sets are thus distinguished based on this criterion : 

- The pulsating variable stars (60 objects) : 
For stars with period P < I day, two samples of 
RR Lyrae are identified : the RRc have R21 < 0.4 
(14 objects) and the RRa6 have R21 > 0.4 (5 ob- 
jects). We a dopt the morphological classification 
proposed by AntoncUo ct al. (1986 ) and classify 
as s-Cepheids the stars that lie in the lower part 
of the R21 — P plane, and as classical Cepheids 
the remaining stars. s-Cepheids pulsate in the first 
overtone and classical Cepheids in the fundamen- 

tal mode (see e.g Beaulieu et al. (1995 ); Beaulicu 

fc Sasselov (1996| )). We use the empirical function 
Rfi^P) = 0.4 - (P/30 days) in order to sepa- 
rate these pulsation modes (see Fig. 7). Among 
the five objects that pass the s-Cepheid cut, only 
three belong to the R21 — P and $21 — P distri- 
butions of galactic s-Cepheids and their phase pa- 
rameter $21 is poorly constrained. Besides most 
of classical Cepheids have amplitudes larger than 
ARe > 0.4"'*s (see e.g [Afonso et al. (1999] )). For 
stars with period P > 1 day, three samples are then 
identified : the classical Cepheids have R21 > i?2i* 
and ARe > 0.4™'^8 (6 objects); the s-Cepheids 
have i?2i < -R21' and 3 rad < $21 < 6 rad (3 
objects); the remaining 32 objects are catalogued 
as miscellaneous pulsating stars. 

- The non-pulsating variable stars (372 objects) : 
The remaining objects have similar amplitudes in 
both passbands. We classify them according to the 
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Fig. 6. Phase diagrams (or light curves for the two lowest panels) for typical variable stars of our catalogue in Re 
magnitude. 



following criteria. Algol systems (130 stars, type 
EA) display well-defined eclipses whose secondary 
one has a depth lower than half the primary one, 
and possibly flat light curve between them. EB 
type objects (35 in total) show a secondary eclipse 
equal to half the primary one. The EW type (11 
objects) is characterised by similar depths of the 
two eclipses. The members of a residual sample of 
196 objects do not look like convincing eclipsing bi- 



naries and are catalogued as miscellaneous variable 
stars. 



As emphasised by Udalski et al. (1999a) a large num- 
ber of variable objects show small amplitude sinusoidal 
variations, such as ellipsoidal binary variable stars. 
A contamination of the sample of pulsating stars by 
eclipsing binaries is thus possible. 

Figures 8 and 9 show the location of the selected vari- 
ables in the colour-magnitude diagrams. Also plotted are 
10,000 stars located in the central part of the two fields 
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Fig. 7. Distinction between classical Cepheids, s-Cepheids 
and miscellaneous pulsating stars in the R21 — P and 
$21 — P planes. The individual uncertainties are reported. 
The cloud of points r epresent the Cepheid s observed in the 
LMC (adapted from Afonso et al. (1999 )). The curve cor- 
responds to the empirical function i?2i'(-F') = 0.4— (P/30 
days) used to distinguish between s and classical Cepheids. 

tm550 and gn450. Most of the Cepheids are much brighter 
than our magnitude threshold (cut 2); this is not so for 
RR Lyrae (see Fig. 9, lower panel). Our catalogue is thus 
not complete for this type of variable stars, as already 
mentionned. 

4. The catalogue 

The catalogue is composed of two tables containing ob- 
jects with periods P smaller or larger than 30 days, respec- 
tively. The identifier of each star is given according to the 
recommendations of the lAU Commission 5 in The Rules 
and Regulations for Nomenclature (see the Annual Index 
of A&A). The general acronym used in the catalogue is 
EROS2 GSA followed by J2000 equatorial coordinates 
in the format JHHMMSSiDDMMSS. The remainder of 
the identifier in parentheses gives some information relat- 
ing to the internal organisation of the EROS database : 
gimnn or tmnnn is the name of the field, followed by 
the CCD number and the location on the image follow- 
ing the EROS II nomenclature. The remaining number 
is the star identifier used in the EROS database. As an 
example, J132630-630945(tm5504ml2359) is the name of 
the 12359th star observed in quarter m of CCD 4 in the 
field tm550. The J2GG0 equatorial coordinates of this star 
are 13:26:30.11, -63:09:45.66. 

The equatorial coordinates (J2000) of individual stars 
have been obtained as follows. First, we have inserted the 



suitable WCS keywords into the header of the EROS II 
reference images using the WCSTOOLS package (Mink 
1999). Whenever possible, the cross-identification of each 



star with previously known objects within a 10" search 
radius has been done using the Simbad and Vizier 
databases available at the CDS, Strasbourg. 
The tables contain the following information : 

1. Identifier 

2. Right ascension a (J2000) 

3. Declination S (J2000) 

4. (Re) mean magnitude in EROS-red 

5. ARe amplitude peak to peak in Re 

6. (Ve) mean magnitude in EROS-visible 

7. AVe amplitude peak to peak in Ve 

8. Period in days. Note that periods longer than 30 days 
are given with less accuracy since the time span of the 
measurements does not allow a precise determination. 
Measured periods which are longer than 60 days (i.e 
2 /3 of the observation period) are flagged by writing 
"P > 60 d" and have no warranty to be true periodic 
variable stars. The peak to peak amplitude of these 
stars could be meaningless and the mean magnitude is 
determined with low accuracy. 

For Cepheids and RR Lyrae the results of the Fourier fit 
are given : 

9. Fourier coefficient ratio i?2i; 

10. Fourier coefficient ratio R31; 

11. Phase difference (j>2i (in rad); 

Also given when possible : 

13. Type of variability (C=classical Cepheids, S=s- 
Ccpheids, puls.= miscellaneous pulsating stars, 
EA,EB,EW=eclipsing binaries, misc=miscellaneous 
variable stars, SR=Semi-Regular variables, M=Miras, 
LPV=long period variables); 

14. Name of cross-identified object(s) within a search ra- 
dius of 10". 

The catalogue is planned to be installed at the CDS (see 
also our Web site ittp : //eros . in2p3 . f r/ ). 

A by-product of such a survey is the possibility to 
update the coordinates given in older catalogues. As 
an example, it was found that the well known classical 
Cepheid 00 Cen with the SIMBAD identifier V*00 Cen 
was 20" away from an EROS object, namely J132630- 
630945(tm5504ml2359), as shown in the finding chart of 
Fig. 10. The period and coordinates we report are compat- 
ible with the ones given by the General Cat alogue of 
Variable Stars (GCVS)0 ( [Kholopov 1985| ) (see Table 
^). It seems that an error occurred when the coordinates 
of this particular star were filled in the SIMBAD. The 
light curve of this object is shown in Fig 6. 

We have performed several cross-identifications of 
our catalogue with those previously available, namely 



Queries on this catalogue are possible 



at 



http://www.sai.msu.su/cgi-bin/wdb-p95/gcvs/stars/form 
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Table 4. Period and coordinates obtained on 00 Cen by 
SIMBAD, GOVS and this study. 



Source 
object 
P (days) 
a (B1950) 
5 (B1950) 



GOVS This study 



vi^as not operational at the time of the observations. 
Therefore the comparison can only be carried out on 0.7 
square degree. Furthermore, as pointed above, this com- 
parison is restricted to stars with magnitude Ic < 14. A 



Tr=i= ^ „ T^nlnnr, nnr^r^Ar/^ — ^ ^ „ , to,tial of 1 18 EROS obiBcts, 6 GCVS aud 23 CKS91 ob- 

V* 00 Cen 00 Cen J132630-630945(tni5504iiil235m . . rr., , • 

^2 ggQ5 12 894 jects lie m this region. Ihree objects are common to the 

13-23-06 85 13-23-09 13-23-08 8 GCVS and CKS91 catalogues one. The overlap between 

-62:54:24.9 -62:54.0 -62:54:11 the three catalogues represents only 19 objects. Among 

them one finds the most interesting ones, such as Cepheids 
OO Cen and V881 Cen (see Table|). Some objects show 
a large difference between the two magnitude determina- 
tions. These are long period variables for which the mean 
magnitude is measured on only a part of the whole period, 
and thus ill determined in both surveys. Only 7 known 
variable stars are not recovered by our analysis (see Table 
Conversely, 99 objects of our catalogue are not listed 
by CKS91. All of them are labelled Long Period Variable 
stars (LPV). 

Table 6. Previously known variable stars not recovered 
by this survey. For each of them we give (where possible) 
their type of variability, their period and the reason for 
their absence in the catalogue. 




Fig. 10. Finding chart of the star 00 Cen - J132630- 
630945(tm5504ml2359). The position given by the 
SIMBAD database is shown by the square which is 10" 
wide, while the EROS position is shown by the circle. 



Object 


Type 


P(days) 


catalogue 


comment 


HQ Nor 


EB 


90.9 


GCVS 


too long period 


HY Nor 


Mira 


236 


GCVS 


close to 










the CCDs gap 


UW Nor 


EA 


8.4860 


GCVS 


fails cut 6 


13214-6256 


LPV 




CKS91 


fails cut 6 


13218-6254 


LPV 




CKS91 


fails cut 6 


13248-6249 


LPV 




CKS91 


fails cut 6 


13232-6249 


LPV 




CKS91 


fails cut 6 



the IRAS Point Source Catalog ( [Beichman et al. 
1998|)7|the MSX5 C Infrared Astrometric Catalog 
( |Egan et al. 1996[), the Two Mic ron All Sky Survey 
(2MASS) (ISkrutskic et al. 1997| ), the CKS91 catalogue 
( [Caldwell et al. 1991 ) and t he General Ca talogue of 
Variable Stars (GCVS) ( [Kholopov 1985| ) . A total of 38 
IRAS sources and 220 MSX5C objects have been thus re- 
trieved. The overlap with the available 2MASS catalogue 
exists only for the EROS field gn459, representing 255 
stars. A total of 233 2MASS objects have been thus re- 
trieved among them 37 objects classified as Semi-Regular 
variables. 



et 



An overlap exists with the CKS91 catalogue (Caldwell 
al. 1991). These authors have searched for bright 
Cepheids and other variable stars with / < 14, to- 
wards Crux and Centaurus, during 42 days with less than 
10 measurements per star. A small overlap exists be- 
tween this survey and the GSA fields tm550 and tm551. 
Unfortunately this overlap involves our CCD #2 which 
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Fig. 8. Colour-magnitude diagrams {Ic vs Vj-Ic) for the Miras (represented by dots •), Semi-Regular variable stars 
(diamonds o) and eclipsing binaries (open circles o). The dotted line corresponds to cut 2 on the luminosity of the 
stars. 
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Fig. 9. Colour-magnitude diagrams {Ic vs Vj-Ic) for the classical Cepheids (represented by stars s-Cepheids (dots 
• ), and RR Lyrae (open circles o). 
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Table 5. Previously known variable stars recovered by this survey. For each star we give its name, the original 
catalogue, the type of variability, the period and the I magnitude given (if any) in the catalogue and in this survey. 



Object 


Catalogue 


Type 


P(days) 


I 


EROS ID 


P(days) 


Ic 


OO Cen 


GOVS 


C 


12.8805 


9.96 


J132630-630945 

(tm5504ml2359) 


12.8941 


10.20 


V881 Cen 


GOVS 


C 


- 


10.57 


J132721-630110 
(tm550313554) 


15.2278 


10.97 


V608 Cen 


GCVS 


EB 


1.6287 


12.95 


J132948-630634 
(tm5505ml3922) 


1.7601 


12.05 


CKS91 


CKS91 


EA 


- 


13.29 


J132758-631449 


5.4918 


13.03 


13246-6259 










(tm550519598) 






CKS91 


CKS91 


E 


- 


13.05 


J133311-630023 


0.2141 


13.20 


13297-6244 








1 


(tm5511m5801) 






CKS91 


GKS91 


SR 


- 


12.35 


J132400-632057 


45.83 


12.34 


13206-6305 










(tm550411003) 






CKS91 


CKS91 


LPV 


- 


11.73 


J132427-631542 


P > 60 


11.75 


13211-6300 










(tm5504k5591) 






CKS91 


GKS91 


SR 


- 


10.51 


J132444-632224 


55.00 


10.79 


13214-6306 










(tm550417786) 






CKS91 


GKS91 


SR 


- 


12.85 


J132447-631201 


52.39 


12.77 


13214-6256 










(tm5504k8841) 






CKS91 


GKS91 


LPV 


- 


10.43 


J132454-631437 


P > 60 


10.23 


13215-625 










(tm5504k9938) 






CKS91 


CKS91 


LPV 


- 


11.81 


J132519-631720 


P > 60 


13.00 


13219-6301 










(tm5504113138) 






CKS91 


GKS91 


LPV 


- 


12.51 


J132433-632201 


58.46 


11.96 


13229-6252 










(tm5504ml0330) 






CKS91 


CKS91 


LPV 


- 


12.48 


J132625-630645 


P > 60 


12.41 


13230-6251 










(tm5504nill603) 






CKS91 


CKS91 


LPV 


- 


12.63 


J132714-634338 


P > 60 


10.78 


13238-6227 










(tm550711058) 






CKS91 


GKS91 


LPV 




12.21 


J132719-633644 


P > 60 


11.91 


13238-6254 










(tm550711752) 






CKS91 


CKS91 


M 




12.46 


J132803-632200 


P > 60 


13.99 


13246-6306 










(tm5505110106) 






CKS91 


GKS91 


SR 




12.66 


J132813-630013 


44.80 


12.72 


13247-6244 










J132813-630013 






CKS91 


CKS91 


LPV 




12.37 


J132948-631227 


P > 60 


12.56 


13264-6256 










(tm5510n4312) 






CKS91 


CKS91 


M 




12.65 


J133222-632339 


P > 60 


12.65 


13289-6308 










(tm5513kl2129) 
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5. Discussion 

The motivation of our search was to improve our knowl- 
edge of the distance distribution and of the extinction of 
the microlensing source stars used in papers I and II. In 
the following we use the RR Lyras which are well-known 
distance indicators and have been observed in all six direc- 
tions that we investigated. The GSA fields having a high 
non-uniform absorption, we give only an average redden- 
ing towards our fields, based on EROS data alone. 

For each selected RR Lyrae we estimate the extinc- 
tion A{V) = 2.4 X E(V — I) using the standard e xtinc- 



tion coefficients (|Schlegel et al. 1998| ; |Stanek 1996|) . The 
colour excess E{V — I) is derived from the colour that we 
measured and their intrinsic colour : {V — /)q ° — 0.4 
and {V -I)^'^" = 0.2 (|Alcok et al. 1998|) . The uncer- 
tainty on individual extinctions is estimated to be O.S™'*^. 
This error includes the 0.2™'*s uncertainty on the EROS 
colour measurement (see Eq. (P) and a O.l™'^^ rnagni- 
tude uncertainty on the intrinsic colour. The absolute 
magnitude of RR Lyras is ~ 0.7 1 with a precision 



of 0.07™''g ( pould and Popowski 1998| ). Finally, we esti- 
mate the distance to each star simply by using the relation 
m — M — by, (log(c?) — l) + A. The typical uncertainty on 
the distance is 20%. The left panels of Fig. 11 show the 
obtained extinction A{V) versus the calculated distance. 
The mean extinction of the RR Lyrffi is A{V) ~ 3.6 ± 0.2 
magnitude towards Mus and 3.3 ± 0.1 towards 7 Nor, 
with a dispersion of 1 magnitude. The mean distance of 
the RR Lyra? is ~4.7± 0.3 kpc towards 9 Mus and 5.3±0.2 
kpc towards 7 Nor. The dispersion of the values is 1.4 kpc 
which reflects the spread in distance of disc stars. The disc 
population contributes to the sources of the microlensing 
events that we observe in the Galactic plane. A model 
based on the distribution of matter in the disc and the 
luminosity function of neighbooring stars has been used 
in Derue (1999 ) to estimate the distance of the star popu- 
lation from the disc. The distance of RR Lyrae is in good 
agreement with the one obtained with this model. 

6. Conclusion 



> 
< 



i5i Mus 



Mean 

RMjS| — I 



4.687 
1.425 



Distonce (kpc) 



Distance (kpc) 



Nor 



+1* 



Mean 

RMS^ 



2 4 6 8 10 

Distance (kpc) 



1.75 
1.5 

1.25 



0.25 - 

L 



5.333 
1.403 



2 4 6 8 10 

Distance (kpc) 



Fig. 11. Extinction A{V) versus the distance of RR 
Lyras (in kpc) towards 9 Mus and 7 Nor. The right panels 
show the inferred distance distribution. 



found it to be ~5 kpc. Yet the statistics being quite lim- 
ited, we are considering pursuing this effort by launching 
a longer search for variable stars based on the results of 
this first campaign. 
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Harvard. The Skycat/Gaia tool is the result of a joint effort 
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